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C^ I interactions of vector mesons and baryons. We estimate the parameters A, /x, u that characterize 



the angular distribution of the produced dileptons in a region of q^ << Q^, where perturbative 
QCD presents an effective strong coupling due to large logarithm corrections. 



X 

^ . PACS numbers: 11.25.Tq ; 13.85. Qk; ll.25.Mj; 12.40.Vv 



* Electronic address: 'c.a.ni.ballonbay ona@durham.ac.uk| 
^Electronic address: ,boschi@if.ufrj .br 



■f Electronic address: braga@if.ufrj.br| 



I. INTRODUCTION 



Proton-proton collisions provide a very large amount of information important for the 
understanding of fundamental interactions. These collisions can produce a large variety of 
particles through many different processes involving strong and electroweak interactions. 
Among the observed particles, one can easily distinguish lepton pairs (a lepton and an 
anti-lepton) usually called dileptons, that can be produced in different ways. Drell and 
Yan proposed a mechanism, within the parton model, to explain production of dileptons 
in hadronic collisions [l|, |2| . This mechanism consists on the annihilation of a quark from 
one hadron with an anti-quark from the other. This annihilation leads to the production a 
virtual particle (typically a photon) that decays into the lepton pair. 

The analysis of dilepton production is important to understand the internal structure of 
hadrons. From the corresponding cross sections one can obtain structure functions related 
to the distribution of partons inside the hadron. 

The differential cross section for dilepton production has an angular dependence that 
can be characterized by three scalar parameters: A, fi, v. These angular parameters have 
been studied using perturbative QCD and the parton model originally in |3H9|. For a more 



recent discussion see 



lOl lll| and references therein. These parameters 



lave been measured 



recently for the case of dimuons by FNAL E866/NuSea Collaboration 12j using data from 
the collision of 800 GeV beams of protons against a hydrogen target. Different hadronic 
collisions lead to different angular parameters. See for instance the case of proton-deuteron 
in pj. 

The perturbative QCD calculations for the Drell Yan cross section work well when the 
transverse momentum of the dilepton is large with respect to the dilepton mass. However, 
when the transverse momentum is small, the perturbative series involves the product of the 
(small) coupling constant as with large logarithm corrections. This combination leads to a 
series with a strong effective coupling. So, the standard perturbative approach is not reliable 
in this regime. In this case one needs a resummation to all orders, as discussed in 11|. 

Recently, alternative approaches to gauge theories at strong coupling were developed 
based on gauge string dualities inspired in the AdS/CFT correspondence |14l4l6l|. In partic- 



ular these dualities lead to some holographic mode 
QCD known as AdS/QCD (see for instance 
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s to describe non perturbative aspects of 



24,l26|,l27|). 



In this article we propose a holographic model to calculate contributions to dilepton 
production through the decay of a virtual photon in proton-proton collisions. Inspired by 
vector meson dominance, in our model the virtual photon giving rise to the dilepton comes 
from the decay of a vector meson. This vector meson is produced by the annihilation of 
two other vector mesons emitted by the protons. The dynamics and interactions of baryons 



and vector mesons are described using the AdS/QCD hard wall model |17H20|. Hadrons 
correspond to modes of a Kaluza-Klein expansion of five dimensional fields living in an AdS 
slice. The size of the slice represents a mass gap in the 4-d effective theory. The hadronic 
masses are determined by the five dimensional wave functions and the boundary conditions 
while the effective coupling constants arise from the integration of interaction terms in the 5- 
d action. These masses and couplings are used to calculate the scattering amplitude relevant 
for the process of dilepton production. For simplicity, we consider only final hadronic states 
with spin 1/2. Using this model we find the angular distribution parameters A, /i, z/, for 
kinematical regimes where the dilepton transverse momenta are small. 

In section 2 we review dilepton production in proton-proton collisions. In section 3 we 
calculate, within the hard wall model, the fermion and vector meson masses and couplings 
relevant for our model. We present in section 4 our model for inclusive dilepton produc- 
tion and estimate the parameters A, fi, z/, that characterize the angular dependence of the 
differential cross section, for kinematical regimes compatible with those analysed by FNAL 
E866/NuSea Collaboration. 

II. INCLUSIVE DILEPTON PRODUCTION IN P-P COLLISIONS 

The production of dilepton from a proton-proton collision through the decay of a virtual 
photon is represented in Figure 1. Two protons with initial momenta Pi and P2 interact 
producing a time-like virtual photon with momentum q plus some additional hadronic states 
which are not observed. The virtual photon decays into a lepton and anti-lepton with 
momenta ki and fco. This section is based on the extensive discussion of dilepton production 
presented in ref . |8| . For a more recent study of angular distribution in Drell Yan process, 



see also 
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The differential cross section, for the unpolarized case, can be written as 



^ d^ki d^k2 



^^ = J^s ^''^-^ 



(27r)32|A;i| {2Tif2\k2 



(1) 



where s = —{Pi + P2Y and g^ = —m'i with m^ > the virtual photon mass. The lepton 
masses were neglected with respect to their momenta and the proton masses were neglected 
with respect to the center of mass energy ^i. 

The hadronic tensor W'^'^ is expressed in terms of the matrix elements of the electromag- 
netic hadronic current as 29 1 



^"^ = \J2J2J d'xe-'^-^Pu Shi, P2, Sh2\ [J^ix), J^x)] |Pi, Shi, P2, Sh2) , (2) 

Shi Sh2 

where Shi and Sh2 are the spins of the initial hadrons. On the other hand, the leptonic 
tensor L^i, is defined in terms of the leptonic current as 
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where Sli and Sl2 are the spins of the leptons and rj^i, = diag(— , +, +, +). 

The most general tensor W^" that satisfies hermiticity, parity constraints and gauge 
invariance can be decomposed as 
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(4) 



where Wi, W2, W3, W4 are the invariant hadronic structure functions that contain the rele- 
vant information for the dilepton cross section. These structure functions depend on four 
scalar variables that are combinations of the momenta Pi, P2 and q. 
It is useful to introduce the helicity structure functions 

1 



Wt = lVi,i , Wl = Wo,o , Wa = ^(^0,1 + ^^1,0) , Waa = W^i,-i 



(5) 




FIG. 1: Illustrative diagram for a Drell-Yan scattering. 



where 

w.,., = vlX'W^'^ ^' ^' = (-^' °' ^) ' (^) 

and r]^ (g) are the polarization vectors of the virtual photon, in its rest frame, defined in 
terms of Cartesian unit spatial vectors X^ = (0, X), y = (0, F), Z^^ = (0, Z) as 

1 



'(0) ' '(±1) a/2 



(7) 



The hadronic tensor can also be decomposed in terms of the helicity structure functions 



as 



W^"" = ir] 



HU 



q^^q" 



{Wt + W^aa) - 2X^X^W^AA 



+ Z^'Z'^iWL -Wt- PVaa) - {X^Z'' + X''Z^')W/^ 
so that the differential cross takes the form (in the photon rest frame) 

^4 



(8) 



da 



]{2Tifm'is 



Wt{1 + cos^ 9) + Wl{1 - cos^ 9) + W^ sin 29 cos , 
+ VFaa siii^ ^ cos 20 sin 9d9d(f) d'^q , 



(9) 



where ^ and are the spherical angular coordinates for the vector ki with respect to the 
Cartesian system X,Y, Z : ki = \ki\ (sin 9 cos 0, sin 9 sin 0, cos 9) . 

In order to compare with experimental results it is interesting to introduce parameters 
that characterize the angular dependence of the cross section. These parameters are defined 
by 



Wt - Wl Wa 

A = TTT TTT- , /i = TT^ TTT- , l^ 
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so that 



v 



rfa ~ 1 + A cos 9 + jj, sin 29 cos + — sin 9 cos 20 . 



(10) 

(11) 



Note that the hehcity structure functions depend on the choice of the unit vectors X, Y, Z 
which must be defined in terms of the momenta Pi, P2 and q. There are different possible 
choices as discussed in f8'|. Here we will follow the CoUins-Soper frame |9|, which is 



+ 02U/2\ 



^/s{m'? + q^) 



Qp 



pi^_ 



Pi-q 



+ qp 



Pf 



Pi-q 



P^ 



q' 

P2-q 



q" + P^ 



P2-q 



q' 



XM 



m^ 



qp 



pf 



-q' 



Pi-q 



+ 



^/sqT{m'^ + qxY^'^ ^" L " q^ 

Pi - ^^^^ - P^^ + ^-^-^"^ 



q- 



P2-q 



q^ 



-q' 



-'2 



P2-q 



Y^" = (0,f) = {0,Z xX) 



(12) 



where 



qp 



q-{Pi-P2) 



qp 



q-{Pi + P2 



qT={ql-q., 



p "^7 



2N1/2 



(13) 



Using the coordinate choice given by (TT2|) in eq. ([8]) and comparing with the expansion 
of the hadronic tensor in terms of the invariant structure functions given in eq. (|1]) one finds 

g2 qlW2 + qpqpW:i + qlWi 



m2 +g|, 
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(14) 



Using perturbative QCD and the parton model, one can calculate the angular distribution 
parameters A and u for large qT, as functions of q and qr, finding in the Collins-Soper frame 
(in our signature where g^ < 0) |5| 



A 



_ g' + jgT 

Pe^* „2 3„2 ' 

q 2^7 



Z/, 



9t 



pert 



-q^ + hr ' 



(15) 



The angular parameter fipert can not be written simply as a function of q and qx, even for 
large qx- Rather, it involves integrals over parton distribution functions, which are not a 
priori known 11 1. 
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In the following sections we develop a holographic model for calculating the invariant 
structure functions Wi,W2,Wz,Wi and then we estimate the parameters A,/i, z/ for kine- 
matical regimes with small qt- 

The kinematical regimes that we will investigate are in the region analyzed recently in 



ref. [12|], considering dimuons produced in collisions of 800 GeV beams of protons against 
a hydrogen target. In this reference de range of dimuon masses analysed was 4.5 < m^^ < 
15GeV and the transverse momenta is in the region < g^ < 4 GeV and < xp < 0.8. The 
mean values found for the angular distribution parameters were: 

(A) =0.85 ; (/i) = -0.026 ; (z/) = 0.04. (16) 

In order to compare our results with those of ref. [12l| we identify m^ = m^^. 

III. VECTOR MESONS AND BARYONS IN THE HARD WALL MODEL 

In this paper we consider the production of dileptons in proton-proton collisions through 
the exchange of vector mesons. So we need first to describe vector mesons and baryons 
and their interactions in the hard wall model. For simplicity, we consider that each final 
hadronic state is just a spin 1/2 baryon. For previous discussions of hadrons in AdS/QCD 



models see, for instance. 



21 



24|. 



The hard wall model consists on a 5-d theory living in an AdSs slice with metric 

ds^ = ^-dt^ + dx^ + dz^] , (17) 

z^ 

where Q<z<zo = l/K and A is an IR energy scale for the dual 4-d effective theory. The 
physical spectrum of the hadronic particles is obtained after imposing boundary conditions 

aX z = zq. 



A. Wave functions of vector mesons 

Consider the action for the non-Abelian 5-d gauge fields in the presence of a gauge-fixing 

term 

r { ^ ....1/1 .. \^1 

(18) 




where Fmn = <9Af^Af— f^Af^j\/+[^A/, ^a^], with Am = A^jT"-, and T°- are the generators of the 
SU{Nf) flavour group. The constant n will be related in the following to the normalization 
of the gauge fields A^^. As we are going to show at the end of section IV, this constant will 
not contribute to the angular parameters A, fi, and z/. 
For the AdS metric (TT7|) and setting A^ = we find 



S = K d!^xdz\J—gh'^ Tr 



^Tl^-V^'F^^F^p - ^rjf^'^d^A^d^A, - hri^-^d^A,) 



(19) 



where h = h{z) = z^ /R^ and rj^^ is the 4-d Minkowski metric. 

Since the z-coordinate is compact, we consider a Kaluza-Klein expansion for the gauge 
field 

oo 

A,{z, x) = f{z) a,{x) + ^ riz) v;{x) , (20) 



n=l 



where the modes /", f^ satisfy Neumann boundary conditions at the IR cut oS z = zo and 
regularity conditions /^(-z) = 1, f"'{z) = at the spacetime boundary z = 0. 

It is important to remark that the infrared Neumann boundary condition arises from the 
gauge invariant boundary condition F^fj, = 0. The latter is a necessary condition to preserve 
the U{Nf) gauge symmetry of the vector field in the bulk. An infrared Dirichlet condition 
is not allowed because it breaks the bulk gauge invariance. For a discussion see ref. 251 ]. 

Imposing the conditions 
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K / dz^^h'r{z)r{z) = 5' 

'o 



1 



-gh^ 



94v^/i29,r(^)] = -minz) , 



(21) 
(22) 



for the modes n = 1, 2, 3, . . ., and the condition 



d,[^^h^dj\z)\=Q 



(23) 



for the non-normalizable zero mode, we obtain a 4-d effective action 
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(24) 



where Smt. represents the interactions terms, we defined 



2/ Y-0N2 



do = K dzy/—g h (/' 



2 /?0 rn 



dn = n i dz^/^h ff 



(25) 



and a^j, = d^a^, — d^a^^, v^^^ = (9^f" — d^v^ are the kinetic parts of the 4-d gauge field 
strengths. 

In order to diagonahze the kinetic terms we redefine the vector field 



K{x) = K{x) - dna^{x) 



(26) 



so that the action takes the form 



/oo , 



^^-^{8^ -ml)+[--l] T^^'^ri'^Pd^d, 



v]^ + gv-v^^vyu 



~^Jint. + 



(27) 



with g^n = m\dn and the dots represent the divergent terms arising from a^. The fields 
■u" are interpreted as vector mesons with mass m„ while the field a^ is interpreted as the 
photon which decompose into vector mesons with decay constant gr^^. This way the hard 
wall model realizes vector meson dominance. 

Note that the operator between brackets in the vector meson kinetic term depends on 
the parameter ^. The 4-d vector meson propagator is given by the inverse of this operator 
in momentum space and takes the form 



A^^{P,ml) 



P"^ + nil 



n 



/ii/ 
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P2 + ^ml 



(28) 



In order to calculate the masses and couplings appearing in eq. (I27jl . we now consider 
the solutions of eq. (!22|) with the chosen boundary conditions 



/"(z) = CnZJi{mnZ) 



(29) 



where rrinZQ are zeros of the Bessel function Jq{w), implied by the Neumann boundary 
condition over f"'{z) aX z = zq. The normalization condition ( 1211) implies that 



(30) 



kR ZQ\Ji{mnZQ)\ 

The solution for the zero mode eq. fl23|) satisfying the corresponding boundary conditions 
is f{z) = 1. 



The coupling constant (^^n takes the form 

Jo 

= KRMm- d,r{z)\^_^ = KRCnTUn. (31) 

Now we consider the interaction Lagrangian for three vector mesons 
Svvv = -K fd^xdzy^h^r]^'''r]''^Ti{[A^,A,]daA(,} 

= Yl ^.".-.^ /"c?^a;r/'^V^Tr {[{};,<] d^v'^} + . . . , (32) 

n,m,i 

where the triple coupling is 

9v-.^.^ = K j dz^^gh{zfr{z)r{z)f{z) , (33) 

and the dot terms are divergent terms arising from the photon a^. 

In order to calculate the masses and couplings for the vector mesons, it is important to 
discuss their dependence on the parameters of the model. The masses depend on the IR scale 
A = l/^o (the same will happen in the fermionic case). We will fix this scale A in section 
IV using the mass of the p meson. On the other hand, the couplings carry a dependence on 
the product kR. In particular, using eqs. ( 130|) and (!3T]) . one finds that the couplings g^-n. are 
proportional to yj kRIs?, while from eqs. f l29|) , fl30|) and f l33|) one concludes that the triple 
couplings g^n^m^e contain a factor l/y/nR. The factor kR will appear as a multiplicative 
factor in the hadronic tensor. So, as we will discuss in section IV, it will not contribute to 
the angular parameters, that involve only ratios of structure functions. So, we do not need 
to fix a value for this quantity. We show in the appendix A some numerical results for the 
masses and coupling of the vector mesons, up to these factors. 

It is important to remark that the approach to vector mesons in the hard wall model that 
we considered here is very similar to the one presented in |26|, |27| within the D4-D8 brane 
model. 

B. Wave functions of spin 1/2 baryons 

In order to describe spin 1/2 states within the hard wall model, we start with the 5-d 
Dirac action 

Sf = kf j d^xdz^^i){D -M)ip, (34) 
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where ip is a 5-d spinor with mass M and the covariant derivate is defined as 

D = ^b'd, + ^V"'l,d,]-^7\ (35) 

with 7^ the 4-d Dirac gamma matrices and 7^ the chirahty matrix. The constant kf will 
be related in the following to the normalization of the fermionic fields ip and ip. As we are 
going to show at the end of this section, this constant will not contribute to the angular 
parameters A, /i, z/ or to the hadronic tensor W'^'^. 
We consider the Kaluza-Klein expansion 

00 

^ = $^[0"Wn + 0"(^)^-]«"(x) 

ra=l 

00 

^ = ^«"(x) \r{z)V.+r{z)V+] , (36) 



n=l 



with V± = (1/2)(1 ± 7^). Imposing the conditions 



-20 ^ r^o y 



^FJ y^-r{z)r{z) = KF J ^/3^-0"(^)0™(^) = <5„^, (37) 



[zd, - (2 + MR)] 0" = -zM, 



n'T ) 



[-zd, + (2 - MR)] 0" = -zM„0" , (38) 

we find the 4-d effective action 



00 „ 

n=l ^ 



(39) 



where M„ are the masses of the baryonic states of spin 1/2 in the 4-d theory. Mi is identified 
with the proton mass and M„ with n = 2,3, . . . , correspond to excited states. 
The normalizable solutions to the eqs. (l38l) are 

riz) = N^z'/'jMR-l/2{MnZ), 

riz) = N^z'/'JnR+i/2{Mnz). (40) 

According to the AdS/CFT correspondence, the dimension of the boundary fermionic oper- 
ator is related to the mass of the bulk fermionic field by 

A = MR + 2 . (41) 
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Baryon states in QCD are associated to a fermionic operator with dimension A = 9/2. For 
this reason we fix the 5-d mass to MR = 5/2. 

We have two possible boundary conditions in the wall z = Zq. We can set 0"(-2o) = 
so that the baryon masses are given by M^Zq = Xn where Xn are the zeros of the Bessel 
function Jmr-ii2- From the normalization condition we obtain 



V KpR'^ Zo\JMJi+i/2{MnZo)\ 

Alternatively, we can choose 4''^{zo) = so that the masses are given by M„zo = Xn where 
Xn are the zeros of Jmr+i/2- The normalization constants in that case take the form 



V KpR'^ ZQ\JMR-l/2{MnZo)\ 

In this work we choose the boundary condition (p'^i^zo) = to calculate the structure func- 
tions and the angular parameters A, /i, v for several kinematical regimes. We also estimate 
the error associated with this choice of boundary condition by calculating A,/i, v for some 
particular kinematical regimes with the alternative condition 0"'(2;o) = 0. 

C. Interaction of baryons and vector mesons 

We can describe the interaction of two fermions and one vector meson considering the 
5-d action 

Sffv = kf d^xdz^f^^'q>"'i)^^A^^ (44) 

that comes from imposing invariance of the action ( IMj) with respect to 4-d gauge transfor- 
mations. Using the Kaluza-Klein expansions for the fields we find 

Sffv = E / ^'^ ^''^ ^>/^ [C^u^'^+ + C-.-^^^^-] ^>' + • • • ' (45) 



n,m,i ' 



where 



^« z 



9„^l = I^F / dZy/^-<j/\z)r{z)<P\z) 



R 

z ~ 



^- ,™„. = «:f / dz^^g-<P''{z)r{z)<P\z) . (46) 



R 



From eqs. (HOj) and (H2l) one finds that the above couplings are actually independent of 
the parameter up. Thus, kf does not contribute to any of our results and does not need 
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FIG. 2: Feynman diagram for a Drell-Yan process mediated 
by vector mesons within our model. 

to be fixed. On the other hand, the fermionic couphngs fH6|) are proportional to 1/\/kR. 
However, as we already mentioned in the case of vector mesons, we do not need to fix this 
quantity either, since it will not contribute to the angular parameters A, /x, v. 

We show in appendix A some of the numerical results for the masses and couplings of 
vector mesons and baryons that we used in our calculations. 

IV. HOLOGRAPHIC DESCRIPTION OF DILEPTON PRODUCTION 

We calculate the contribution to dilepton production represented in the Feynman diagram 
of Figure 2. This corresponds to the interaction of two protons with momenta Pi and P2 
through the exchange of vector mesons f" and f™. These vector mesons combine into 
another vector meson v^ that decays into a time-like photon that eventually gives rise to 
a lepton pair. At lowest order, the final state corresponding to each proton is one excited 
baryon. These final baryons are not measured, so the hadronic tensor W^'^ involves the sum 
over all possible baryonic states X and Y . Here we will consider only final states X and Y 
with just one spin 1/2 baryon each. 

A. Scattering amplitude 

The optical theorem relates the total cross section of proton-proton-photon scattering to 
the forward scattering amplitude represented in Figure 3. As a consequence, the structure 
functions W^i, W2, W^a, W4 associated with the hadronic tensor W^^ in eq. (jl]) can be obtained 
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P2 Py P2 

FIG. 3: Feynman diagram for the proton-proton-photon forward scattering 



from the imaginary part of the forward scattering tensor T^'^. 

From the diagram of Figure 3 and the Feynman rules that come from the 4-d effective 
Lagrangians of the previous section we find the ampUtude 
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(47) 



where u\-. and -uL) ^^^ the spinors representing the initial protons with momenta Pi and P2 
and 



A''P{P,ml 



P2 + mf 



v^'-ii-O 



pcrpp 

P2 + ^mf 



Ca'P'T — Va'l3'{Pl — Px — P2 + Py) 
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+ VP'AP2 -Py + q)a' - Va'riPi - Px + q) p> (48) 



Note that the denominators of the fermionic propagators do not appear in eq. fl47j) since 
they are transformed into delta functions when one takes the imaginary part of the forward 
scattering tensor. These delta functions impose on-shell conditions on the momenta Px = 



{Ex, Px) and Py = (Ey, Py), implying that Ex = ^^ P^ + M^^ and Ey = ^ P^ + M^^ 
Summing over the spin of the initial hadron 1 we find 



■Shi 



PlPx ^s"x^,^„l^s"x^,n„l Tr (7A7a7a7a^+) + ^s"x^,.->„l^^z-x„n„l Tr (7A7a7a7a^- 



+ -MiM„^ 



,+ 



where the traces over gamma matrices are 



Tr (7<57a) (49) 



Tr (^7A7a7a7a'^±j = '^[VXaVaa " VXaVaa + VXaVaa] + ■■■ 

Tr (-fala) = -^Vaa (50) 

and the dots represent the imaginary term ^2iexaaa whose contribution is negligible. For 
the sum over the spin of hadron 2 we have a similar result. 

In order to perform the integration in the three momenta Px and Py we choose a frame 
where Pi + P2 — g = 0. That means, we work on the center of momentum frame of fi- 
nal hadrons. However, it should be stressed that the results for the structure functions 
Wi, W2, W3, W4 are frame independent. Then the integral takes the form 



2Ex J lEy \ / \ / 

= / , "^'^-^ - 5(Ex + JeI-MI^+MI^-E^^... (51) 

where Ep = Ei + E2 — qo is the energy available for the final hadrons. We can express the 
three momentum volume in terms of the energy as 



d'Px = Ex J El - M2 sin 9xdExdexdct>x • (52) 
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Here Ox and ipx are spherical angular coordinates for the vector Px- This way, the delta 
function can be integrated leading to the condition 

Ex = ^^{eI + MI^-MI^). (53) 

Then the momentum integrals in eq. ( H7|) become just angular integrations, i.e. 



,, d9x dct>xsm9x^ (•.•)• (54) 

4 Jo Jo Ef 

The integrand depends on the internal momenta Px and Py (in the Feynman graph of 
Figure 3 ). In the center of momentum frame of final hadrons, these momenta are expressed 

as: 

Px = (-E'x, l-Pxl sin^x COS0X, iPxl sin6'x sin0x, |-Px| cos6'j!c) 

Py = (Ef-Ex^-Px). (55) 



where \Px\ = \ E\ — M^ . The integrals in 6x and (px are performed for each term in the 
sum over vector mesons and fermions. 

The amplitude defined in eq. (147|) involves sums over all intermediate vector mesons 
and final baryons. For the final baryons X and Y, there is a physical condition that sets a 
natural cut off for the sums, since they are on shell particles. The available energy for them 
is Ep = El + E2 — qo, so, the final states must satisfy: 

M„^ + Mny < Ef . (56) 

For the vector mesons, there is no physical cut off since they are off-shell. In our numerical 
approach, we used the convergence of the sums as the criterion to define how many vector 
mesons should be summed in the internal lines. We computed the sums over the indices 
m, n, m, n from 1 to some integer value rimax starting with rimax = 8 and increasing its value. 
We did the same for the indices £, i taking the corresponding sums from 1 to some integer 
nmax2- We found a good convergence at Umax = 15 and nmax2 = 25. More precisely, adding 
one more term to each vector meson sum, that means taking rimax = 16 and nmax2 = 26 we 
found a relative variation smaller than 10~^ in expression (H7|) . So, we used Umax = 15 and 
nmax2 = 25 in all our calculations. The relative error in the angular parameters, associated 



with this cut off, is of order 10 ^. 
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Kinematical regimes 


Our results 


Perturb 


. results 


V 


12 


93 


V^ 


qr 


XF 


A 


^ 


V 


^pert 


VpeH 


25.8 


1.61 


11.3 


38.8 


1.47 


0.373 


0.76 


0.069 


-0.084 


0.960 


0.0201 


25.6 


1.29 


10.8 


38.9 


1.20 


0.361 


0.76 


0.066 


-0.074 


0.973 


0.0136 


24.8 


0.937 


9.69 


38.9 


0.900 


0.331 


0.73 


0.032 


0.069 


0.985 


0.00770 


24.1 


0.565 


8.40 


38.9 


0.601 


0.294 


0.70 


0.034 


-0.0024 


0.993 


0.00346 


22.7 


0.162 


6.02 


38.9 


0.302 


0.218 


0.57 


0.019 


0.030 


0.998 


0.000876 



TABLE I: Angular parameters calculated from our model and from the perturbative approach, for 



mz 



104 GeV^. The kinematical variables p,q2,Q3,\^,QT are expressed in GeV. Note that /ipert 



is not known, since it depends on the unkown parton distribution functions jll||. 



In 



In our amplitude ( H7|) we wrote for the propagators of all vector mesons of momentum P 
and mass rrii the expression 



A^-'iP,m^) 



P^ + mf 



v^^-a-o 



ptipv 



P2 + ^rnf 



(57) 



This propagator depends on the 5-d gauge fixing parameter ^ so depending on the choice 
of ^ we may have different effective 4-d models. There are three common gauge choices : 
i) .^ — 7- oo which leads to transverse currents but is hard to implement at high energies, ii) 
^ = 1 which does not lead to transverse currents unless we turn on Goldstone bosons and 
iii) ^ = which leads to transverse currents and has a nice high energy behavior. 

As we explained above, the calculation of the dilepton cross section involves momentum 
integration and several sums over fermion and vector meson masses. For this reason, we 
found convenient to work with the gauge ^ = where the propagator leads to transversality 
of the hadronic tensor and has a good numerical behavior in the sense that it reduces to 
a massless propagator in the limit of large P^. We also made some numerical tests with 
the gauge ^ — t- oo but in that case the mass dependence on the numerator implies a huge 
running time making very hard to guarantee convergence of the different sums. 

B. Numerical set up and results 

We fix the size of the hard wall AdS slice zq = 1/A, by fitting the mass of the p meson 
nip = 0.776 GeV. This gives A = 0.323 GeV . Concerning the parameter kR that, as 
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Kinematical regimes 


Our results 


Perturb 


. results 


V 


92 


93 


^ 


qr 


XF 


A 


M 


u 


■^pert 


VpeH 


25.8 


1.53 


11.3 


38.9 


1.50 


0.398 


0.60 


0.11 


-0.15 


0.925 


0.0377 


25.5 


1.20 


10.8 


38.9 


1.20 


0.386 


0.61 


0.085 


-0.11 


0.951 


0.0246 


24.8 


0.859 


9.69 


38.9 


0.900 


0.355 


0.57 


0.070 


-0.038 


0.972 


0.0141 


24.1 


0.501 


8.40 


38.9 


0.600 


0.316 


0.49 


0.057 


0.025 


0.987 


0.00633 


22.6 


0.0808 


6.02 


38.9 


0.302 


0.236 


0.32 


0.021 


-0.013 


0.997 


0.00162 



TABLE II: Angular parameters calculated from our model and from the perturbative approach, 
for mi = 56.3 GeV^. The kinematical variables p,q2,Q3,\^,QT are expressed in GeV. Note that 



fipert is not known, since it depends on the unknown parton distribution functions 



Hi 



1.2 r 

1-1 : 

1.0 : 

0.9 

^ 0.8 : 

0.7 : 

0.6 

0.5 : 

0.4 L^ 
0.2 



ml = 104 GeV^ 



0.4 0.6 0.8 1.0 1.2 

IT (GeV) 



ml = 104 GeV^ 



0.2 0.4 0.6 08 1.0 1.2 



qr (GeV) 



mi = 104 GeV 



0.00' 

0.2 0.4 0.6 0.8 



qr (GeV) 




FIG. 4: We plot our results (blue lines) and the perturbative predictions (red lines) for the angular 
parameters A,/x, i^ for mi = 104GeV^. Note that /i is not compared with the corresponding 
perturbative result fipert since the later it is not known, thanks to its depence on the unknown 
parton distribution functions [lll |. 
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ml = 56.3 GeV^ 



ml = 56.3 GeV^ 




0.4 0.6 0.8 1.0 1.2 

qr (GeV) 



0.2 



0.00 ^ 
0.2 



0.8 1.0 

qr (GeV) 



: 56.3 GeV^ 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 

qr (GeV) 



FIG. 5: We plot our results (blue lines) and the perturbative predictions (red lines) for the angular 
parameters A, fi, v for m^ = 56.3 GeV . Note that [i is not compared with the corresponding 
perturbative result [Xpert since the later it is not known, thanks to its dependence on the unknown 



parton distribution functions 



Hi 



discussed in section III, shows up in all the couplings, one can see from the form of the 
hadronic tensor eq. WI\ that it will contribute just to a common multiplicative factor. 
So, the same factor will appear in all the structure functions. In particular, the helicity 
structure functions defined in eq. ([8]) will share the same multiplicative factor of kR. Thus, 
the parameters A,/i, z/, defined as ratios of linear combinations of these objects, will be 
independent of kR and we do not need to fix it. 

In our frame, of the center of mass of the final hadrons, the protons and the photon 
momenta can be written as 



P2 



p2 + M2,0,0, p) ; q={^m^^ + ql + ql,{],q2,q^) 
ql + (93 - VY + Mf ,0,q2,q3-p), 



(58) 



where Mi = 3.83A = 1.2 GeV is the "proton" mass in the model, once the infrared scale was 
fixed by the mass of the p meson. We analyzed different kinematical regimes, contained in 
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the region analyzed from experimental data in ref. [12|]. Each one was defined by a choice of 
the values of p, q2, qs and m^. From these quantities we obtained the corresponding values 
for a/s, qt and xp = "^^-^ where qx and g^ were defined in eq. (IT^ . 

The numerical calculation where performed using the package "Mathematica" . From eq. 
( H7|) we computed the diagonal elements of the hadronic tensor W^'^. Using these results it 
is possible to invert eq. (jl]) and obtain the structure functions Wi, W2, W3, W^. Then from 
eq. flT4|) we obtain the helicity structure functions in the Collins-Soper frame from which 
we calculate the angular parameters A,/i, z/ defined in eq. ( TTOj) . 



We chose kinematical regimes in the region analyzed in ref. 12|. They covered the range 



20 < m^^ < 225 GeV^, so we decided to test, as representatives of this region, the values 
w? = Tn? = 56.2 GeV^ and m'i = 104 GeV^. For the transverse momentum qx they 
covered the region < gr < 4GeV, so we chose values in the non perturbative range 
0.3<gT< 1.5 GeV. 

In Tables [l] and [TTl we show our choices of kinematical regimes and the corresponding 
results for the angular parameters from our model and from the perturbative expressions 
for A and u ( IT5l) . since the perturbative expression for fi depends on the parton distribution 
functions which are not known. In Table [l] we fixed mi = lOOOA^ = 104 GeV^ while in 
Table [TTl ^ = SSOA^ = 56.2 GeV^ 

We also present plots comparing our results for A, /i, and u with the perturbative ones 
Xpert and Upert in Figures 4 and 5. We see that our model predicts a decrease in the value of 
A as qT decreases. This contrasts with the perturbative expectations of A — )■ 1 when g-r — ^ 0. 
The mean value < A >= 0.85 obtained from the available experimental data is lower than 
the perturbative and greater than our model results. 

Now we estimate the error associated with the choice of boundary conditions for the 
fermions by calculating the angular parameters using the alternative boundary condition 
4>"'{zq) = 0. For this purpose, we calculated again all the fermionic masses and coupling 
constants for this alternative condition. We compare in table III the angular parameters 
obtained using the original boundary condition with those obtained using the alternative 
boundary conditions. From this table we can have some estimative of the relative errors 
associated with the choice of boundary conditions. We find 

6\ 611 5v , , 

— ~ 0.04 -0.06, —-0.6-0.8, 0.06-1.4. (59) 

A ^ v 
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Kinematical regimes 





riginal B.C. 


Alternative 


B.C. 


P 


q2 qs rn^ 


QT 


V~s 


A 


fJ, V 


Aa l^a 


Va 


25.8 
22.6 


1.61 11.3 104 
0.0808 6.02 56.3 


1.47 
0.302 


38.8 
38.9 


0.76 
0.32 


0.069 -0.084 
0.021 -0.013 


0.73 0.13 
0.34 0.048 


-0.079 
-0.072 



TABLE III: Angular parameters calculated with the original boundary condition (/("(zq) = and 



with the alternative boundary condition 
are expressed in GeV. 



1^0 j 



0. The kinematical variables p, q2,qz-, V^, qr, ^ 



These results show that the parameter A has a low sensitivity to the choice of boundary 
conditions for the fermions. This is the largest angular parameter and the one for which 
the predictions of our model are robust and closer to both experimental and perturbative 
results. The experimental results for fi and u are much smaller and present oscillations. In 
our model these parameters are very sensitive to the choice of fermionic boundary conditions. 
However, the absolute values of our parameters are of the same order of the experimental 
ones. 

V. CONCLUSIONS 



As we mentioned in the introduction, the perturbative QCD approach provides a good 
approximation for the dilepton production cross section when the dilepton transverse mo- 
mentum qt is large compared to the virtual photon mass m^. In this approach the hadronic 
tensor W'^'^ depends, through collinear convolution, on the generalized parton distribution 
functions and partonic cross sections. The dominant processes are the Drell Yan quark an- 
nihilation {qq — )■ 7*), the quark annihilation with the emission of a gluon [qq — )■ '~f*g) and 
the quark gluon scattering qg — )■ 7*g (where 7* is the virtual photon). However, in the limit 
of small dilepton transverse momentum {qrp <^ m^) the perturbative approach suffers from 
large logarithmic corrections ll|. Since the data extracted from recent experiments (like the 
FNAL E866/NuSea Collaboration) include this problematic regime in which the effective 
coupling is large it is worth to explore alternative approaches to the regime of small gr- 
in this article we presented a holographic model where dileptons are produced in proton- 
proton collisions by the exchange of vector mesons. Our approach deals only with dileptons 
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that are produced from time-like virtual photons. The model describes only processes where, 
after the collision, each proton transforms into a single baryon of spin 1/2. It would be 
interesting to extend this study in order to include the production of more general final 
hadronic states, like baryons of spin 3/2. 

We applied our model to estimate the cross section parameters A, /i, v that characterize 
the angular distribution of dileptons produced in a proton-proton collision for small qr- 
We compared our results for A and v with those from perturbative QCD. We found that 
in the region analysed A is smaller than Xpert and increases with qr in contrast with the 
perturbative result. For g^ — )■ the perturbative result is Xpert -^ 1 while our model gives 
decreasing values for A. 

As a final remark, it is important to mention that there are other processes that can 
generate dileptons in a proton-proton collision and may contribute to the parameters X, fi,i>. 
Then it would be interesting to explore possible extensions of the present model to include 
these other processes. 
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VI. APPENDIX: TABLES OF MASSES AND COUPLINGS 

Here we list some masses and couplings for vector mesons and spin 1/2 baryons in the 
hard wall model. 

Vector meson masses are obtained from the zeros of Bessel function Jo{w). The g^" 
coupling constants, calculated using eq. ( 13T|) . represent the decay of a vector t>" into a 
photon. The vector meson triple couplings Qv^v^v^ are calculated from eq. ( I33l) . Some of 
these results are shown in Table IIVI 

The masses for spin 1/2 baryons are calculated from the zeros of the Bessel functions 
Ji{w). The triple couplings among two baryons and one vector meson are calculated from 
eq. ( l46l) . Some of these results are shown in Table |Vl 
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n 


A 


9^" 

v';jflA2 


fiVfiKR9„i,i„„ 


BV6Ki?a„i„2,„ 


BVBKi?9„i„3„„ 


6VfiKR9„2,2„„ 


6VfiKi?S„2„3„„ 


6VBKi?a„3„3„„ 


1 


2.4048 


6.5510 


1.0923 


-0.31784 


0.0035072 


0.74819 


-0.34179 


0.71437 


2 


5.5201 


22.943 


-0.31784 


0.74819 


-0.34179 


0.12962 


0.40807 


0.10405 


3 


8.6537 


45.084 


0.0035072 


-0.34179 


0.71437 


0.40807 


0.10405 


0.13996 


4 


11.791 


71.736 


-0.00039171 


0.0053568 


-0.34723 


-0.37659 


0.36944 


0.076647 


5 


14.931 


102.23 


0.000083373 


-0.00074121 


0.0060115 


0.0087557 


-0.38591 


0.32767 


6 


18.071 


136.13 


-0.000024843 


0.00018336 


-0.00090370 


-0.0015000 


0.010128 


-0.39707 


7 


21.212 


173.13 


9.1480x10"'' 


-0.000061090 


0.00023961 


0.00042793 


-0.0018924 


0.011905 


8 


24.352 


212.98 


-3.9015x10-6 


0.000024521 


-0.000084592 


-0.00015811 


0.00057805 


-0.0024352 


9 


27.493 


255.50 


1.8545x10-'' 


-0.000011201 


0.000035649 


0.000068680 


-0.00022572 


0.00079693 


10 


30.635 


300.51 


-9.5825x10-'^ 


5.6307x10-'' 


-0.000016972 


-0.000033399 


0.00010265 


-0.00032869 


11 


33.776 


347.90 


5.2915x10-'^ 


-3.0479x10-® 


8.8399x10-" 


0.000017666 


-0.000051882 


0.00015630 


12 


36.917 


397.55 


-3.0844x10-'^ 


1.7504x10-" 


-4.9345x10-" 


-9.9764x10-" 


0.000028362 


-0.000081999 


13 


40.058 


449.36 


1.8807x10-'^ 


-1.0552x10-" 


2.9110x10-" 


5.9383x10-" 


-0.000016480 


0.000046268 


14 


43.200 


503.25 


-1.1911x10-'^ 


6.6235x10-'^ 


-1.7968x10-" 


-3.6912x10-" 


0.000010057 


-0.000027626 


15 


46.341 


559.13 


7.7939x10-* 


-4.3029x10-^ 


1.1518x10-" 


2.3795x10-" 


-6.3897x10-" 


0.000017263 



TABLE IV: Some masses and coupling constants of vector mesons in the hard wall model. 



n 


A 


BVK-Rot, „ , 








1 


5.1356 


0.92525 


-0.28634 


0.0063502 


1.1995 


-0.21837 


-0.0010294 


2 


8.4172 


0.22261 


0.55996 


-0.29692 


-0.55390 


0.68808 


-0.18644 


3 


11.620 


-0.49828 


0.11371 


0.52173 


0.056561 


-0.62284 


0.59639 


4 


14.796 


0.17584 


-0.49341 


0.078525 


0.087530 


0.19693 


-0.63706 


5 


17.960 


-0.017595 


0.20755 


-0.49104 


-0.070083 


0.040080 


0.25764 


6 


21.117 


0.0053409 


-0.019083 


0.22278 


0.043309 


-0.058728 


0.014012 


7 


24.270 


-0.0022875 


0.0059035 


-0.019476 


-0.029528 


0.036904 


-0.052222 


8 


27.421 


0.0011692 


-0.0026014 


0.0060597 


0.021530 


-0.025401 


0.033015 


9 


30.569 


-0.00066745 


0.0013668 


-0.0027067 


-0.016447 


0.018679 


-0.022778 


10 


33.717 


0.00041141 


-0.00079982 


0.0014436 


0.013002 


-0.014384 


0.016797 


11 


36.863 


-0.00026846 


0.00050392 


-0.00085719 


-0.010551 


0.011455 


-0.012977 


12 


40.008 


0.00018311 


-0.00033522 


0.00054749 


0.0087422 


-0.0093591 


0.010371 


13 


43.154 


-0.00012940 


0.00023256 


-0.00036882 


-0.0073665 


0.0078027 


-0.0085037 


14 


46.298 


0.000094153 


-0.00016685 


0.00025883 


0.0062949 


-0.0066124 


0.0071147 


15 


49.442 


-0.000070205 


0.00012304 


-0.00018766 


-0.0054432 


0.0056801 


-0.0060500 



TABLE V: Some masses and coupling constants of fermions in the hard wall model with the 
boundary condition ^"'(zo) = 0. 
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